Acoustic transient event reconstruction and sensitivity studies with the South 

Pole Acoustic Test Setup 



Jens Berdermann 3 , for the IceCube collaboration b 

"DESY, D-15735 Zeuthen, Germany 
b http://icecube. wise, edu/ 



Abstract 

The South Pole Acoustic Test Setup (SPATS) consists of four strings instrumented with seven acoustic sensors and 
transmitters each, which are deployed in the upper 500 m of the IceCube holes. Since end of August 2008 SPATS 
is operating in transient mode, where three sensor channels of each string, located at three different depth levels, are 
used for triggered data taking within the 10 to 100 kHz frequency range. This allows to reconstruct the position of the 
source of acoustic signals in the antarctic ice with high precision. Acoustic signals from re-freezing IceCube holes 
are identified. All detected acoustic events seen are associated to sources caused by human activities at the South Pole 
111] . Further, the sensitive volume for v interactions outside the IceCube instrumented area has been determined by 
simulation and a flux limit for high energy neutrinos was derived. 

Keywords: Acoustic neutrino detection, Acoustic transient data, SPATS 
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1. Vertex reconstruction 

The reconstruction of acoustic transient events is 
based on the solution of the idealised global position- 
ing equation system 

(x„ - x ) 2 + (y n - y f + (z n - zof = [v s (t„ - t Q )] 2 , 

where four sensor positions and the signal arrival times 
t n ,x n ,y„,z n are used in a single reconstruction n 
I , ... 4 The event vertex is located at the space time 
point to,Xo,yo,Zo- As an idealisation a sound propaga- 
tion in ice without refraction and a constant velocity of 
Vj = 3850 m/s is used. The assumption of a constant 
speed of sound is only suitable for events below a depth 
of around 200 m and leads to a spread of reconstructed 
events for shallower depths, as one can see from simula- 
tion. Solving the equation system above leads to two ex- 
act real solutions, where one of them is the event vertex 
located at the space time point to,xo,yo, Zo and the other 
turns out to be unphysical. The recently used SPATS 12- 
channel configuration [7] allows statistical predictions 
by use of all possible sensor combinations per acoustic 
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Figure 1: Horizontal distribution of all transient event vertices 
recorded since August 2008. The sources of transient noise are the 
Rodriguez Wells (RW), large caverns melted in the ice for water stor- 
age during IceCube drilling, and the re-freezing IceCube holes. 
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event. In case of a noise hit in a sensor the reconstruc- 
tion algorithm for this combination does not converge or 
the result lies far outside the sensitive SPATS area. To 
improve the reconstruction result a cut at the tails of the 
vertex distribution of all sensor combinations has been 
applied. A horizontal distribution of all reconstructed 
transient events recorded since August 2008 is shown in 

Fig.ru 

2. Events from IceCube holes 

Acoustic events were observed from nearly all Ice- 
Cube holes drilled when transient data taking was ac- 
tive. The results for different holes are similar, therefore 
hole 8 1 is chosen below as an example. Events are ob- 
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Figure 2: Events around the region of hole 81 between 10.12 till 
31.12.2009 in the x-z plane. Visible are events before (firn drill), dur- 
ing and after (re-freezing) enhanced hot water drilling starts at the 
20.12. The red line indicates the bore hole position. 

served within 20 days in the hole region (±20 m with re- 
spect to the center of hole 81) during the periods of firn 
ice drilling (< 50 m depth), bulk ice drilling (50 - 2500 
m depth) and re-freezing. Fig. [2] shows the depth posi- 
tion of the events versus the x-position, in chronologi- 
cal order. A few early events observed at 40 m - 100 m 
depth are connected with noise from the firn drill hole. 
During bulk ice drilling events are found in the same re- 
gion but also at larger depth. Strong sound production 
starts about three days after drilling is finished, due to 
the refreezing process. 

About 30 % of the registered events from this hole are 
concentrated in two spots at 120 m and 250 m depth but 
reaching down to about 600 m . The reason is that the 
hole doesn't re-freeze homogeneously but forms frozen 



45 ice plugs between still water filled regions. The pres- 

46 sure produced in this way may give rise to cracks near 

47 the ice water boundary which appear with sound in the 

48 10-100 kHz frequency region. Relaxation later contin- 

49 ues within "arms" freezing upwards to the hole surface 
so and down to the lower ice plug. 

51 Besides providing information about the re-freezing 

52 process of water-filled IceCube holes one can use the 

53 corresponding acoustic events also to understand the 

54 precision of the vertex localization algorithm. The av- 

55 erage values for the (x,y) position of hole 8 1 are deter- 

56 mined to: 40.1+0.1 m in x and 39.4+0.1 m in y. The 

57 width of the distributions is 2.4 m and 4.6 m respec- 

58 tively to be compared with a hole diameter of about 70 

59 cm. The calculated values deviate from the actual ones 

60 by 1.4 m (in x) and 3.9 m (in y). The possible reason for 

61 this deviation will be discussed in the simulation section 

62 E]below. 

63 3. Noise from Rodriguez-wells 

64 When the first acoustic events had been reconstructed 

65 during the period from August to November 2008 (quiet 
ee period) a strong clustering in a certain region of the x- 
67 y plane at about (-150 m, 300 m) became visible. It 
es was found that this was the position of the 2007/08 Ro- 

69 driguez well (Rod-well or RW for short) used for the hot 

70 water drilling system. 

71 This type of well has been introduced by Rodriguez 

72 and others in the early 1960s [2J] for water supply at a 

73 glacier in Greenland. Hot water cycled by a pump sys- 

74 tem is used to melt ice below the firn layer at 60-80 m 

75 depth to maintain a fresh water reservoir. An expanding 

76 cavity is formed with a diameter as large as 15-20 m. 

77 For IceCube and its predecessor AMANDA this tech- 

78 nique has been used in connection with drilling at the 

79 South Pole since mid 1990s. If the well is used a second 
so time a year later, a second cavern is formed at a deeper 
si level. 

82 Having identified acoustic events arising from the 

83 2007/08 Rod-well, three other event clusters were 

84 found, two of them could be attributed to other Ice- 

85 Cube Rod-wells, from 2006/07 and 2004/05-2005/06. 
se The fourth event cluster turned out to be located at the 
87 probable position of the last AMANDA Rod-well used 
ss in the final two drilling seasons up to 2001. No spe- 

89 cific coordinates could, however, be found documented 

90 for that position anymore. The acoustic events from the 

91 Rod-wells used during two seasons are located at larger 

92 depths than those from Rod-wells used only once. The 

93 latter were seen to emit acoustic signals from smaller 

94 and smaller regions around the well core and finally 



stopped, the older one in October 2008, the younger one 
in May 2009. 

In contrast to that, acoustic events are observed un- 
til today from the six and ten years old deeper wells 
(see Fig. [3}. The mechanism of sound production in 
and around the Rod-well caverns is still under debate in 
particular for the older wells. 



Table 1: Parameters of acoustic events from different Rod- 
wells. 
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Figure 3: Acoustic events distribution for y-coordinate versus time. 123 

Lines: fitted positions of AMANDA-RW (red), 04/05-05/06 IC-RW i 24 

(green), 06/07 IC-RW (blue), 07/08 IC-RW (magenta). The zero at 12g 
the time axis indicates Jan 1, 2009. 



02 4. Acoustic event simulation 

103 A simple approach to acoustic transient event 

104 simulation can be done by calculating the sig- 

105 nal propagation times for the distances d n = 
ice -\J(x n - x) 2 + (y n - y) 2 + (z n - z) 2 between source (e.g. 

107 IceCube hole at x,y,z) and sensors n = 1, ... ? n max with 

108 At n = d n /v s . The signal is randomly transmitted from 

109 a certain cylindrical volume (radius 2 m, depth 2000 
no m) around the source. Although knowing that the true 
m IceCube hole diameter is about 70 cm, we take into 

112 account the possibility that tension cracks might ap- 

113 pear outside the hole bounding surface which suggests 
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Figure 4: Reconstructed simulated events for visible holes of the drill 
season 08/09 



a larger simulation radius. The reconstruction of events 
simulated with constant speed of sound and without 
considering attenuation effects implies an exact source 
localization, which is in contradiction to the real data 
vertex results, where a specific data spread around the 
source (FigJTJ and a lack of events below and above 
a certain depth (Figf2j) is visible. The major reason 
for misreconstruction of events at shallow depth (-200 
m < z < m) is probably the depth dependence of 
sound speed [J] which is, therefore, included in the 
simulation. Above 170 m depth, the parametrisation 
v s = -(262.379+ 199.833 |z|i - 1213.08 kl^ms" 1 and 
below a constant sound speed of v s = 3850 ms _1 is used. 
Further improvement is achieved using additional infor- 
mation on sound pressure wave attenuation in the ice 
H. We apply S" R = (S ■ d /d n )e- d " /A with the initial 
amplitude So at the distance do chosen to fit the real 
data (5J and an attenuation length A of 300 m as mea- 
sured for South Pole ice with SPATS 01 • If the signal 
strength at a sensor is above ~ 300 mV (5.2 cr above the 
noise level) a hit is triggered as in real data. 

A good agreement between reconstructed real and re- 
constructed simulated events is finally obtained as one 
can see by comparison of Fig|4]with FigQ]and Figf2] 
As expected we observe a large influence by the depth 
dependent sound speed on reconstructions in the upper 
region of SPATS (between and 200 m depth). The 
significant deviation of the sound speed from the con- 
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142 stant value used in the reconstruction explains the ver- 

143 tex spreading seen in the real data, whereas the direction 

144 of these "smearings" is caused by the detector geometry 

145 and points towards the center of SPATS. The inclusion 

146 of the attenuation length makes it more difficult to ob- 

147 serve deep events which is in agreement with the real 

148 data distribution. 

149 5. Noise from regions outside IceCube 

150 In order to determine the number of events not con- 

151 nected to IceCube construction activities in the sensitive 

152 region of SPATS we omit the area of IceCube strings 

153 and the data-files from the drill periods keeping in mind 

154 that we have a lot of acoustic hits here due to detector 

155 construction. Furthermore we look at depths between 

156 200 and 1000 m, in the region of constant speed of 

157 sound, to avoid the smearing described in section|4] In 

158 the 245 days of transient data taking we found no events 

159 for the recent SPATS- 12 detector configuration in the re- 

160 gion defined above which allows to determine a limit on 

161 the cosmogenic neutrino flux. 

162 In order to determine a flux limit an effective vol- 

163 ume for the SPATS-12 configuration is calculated. Neu- 

164 trino interaction vertices were simulated in a cylindrical 

165 volume with an energy dependent sensitivity radius be- 

166 tween 0.4 - 2 km around the centre of SPATS and at 

167 depth between 200 and 1000 m. We omit again the area 

168 of IceCube strings. The neutrinos were assumed to be 

169 down going with random 9 and <p. Together with the in- 

170 teraction vertex, the direction (9, <p) defines the plane of 

171 acoustic pressure wave. The sensor observation angle 

172 was then calculated relative to this plane for each ver- 

173 tex. 10 7 events were simulated for neutrino energies 

174 E Y increasing from 10 18 eV to 10 22 eV. The hadronic 

175 component of the neutrino energy E^d was assumed 

176 to be a constant fraction y of the neutrino energy i.e. 

177 Ehad = 0.2E V . In this simulation each acoustic sensor is 

178 simply a point at which to calculate the acoustic pres- 

179 sure Pmax as a function of 9 and R. To get a reasonable 
iao value for the signal strength at the sensors we apply the 
iai Askaryan |6[] model assuming a cylindrical energy de- 

182 position in the medium of length L and diameter d. No 

183 attempt has been made to model angular sensitivity or 

184 frequency response. A minimum threshold of ~ 300 

185 mV, like in the real SPATS-12 measurement is used, 

186 which transforms to a necessary minimum pressure of 

187 ~ 70 mPa [7]. The observation of zero events inside 

188 the effective volume of SPATS-12 gives an upper limit 
iao of N bs = 2.44 events by a poissonian 90 % confidence 

190 level il]. The flux limit is shown in Fig |5] together with 

191 the cosmogenic neutrino flux predictions and results of 




log io (E v )[GeV] 



Figure 5: Neutrino flux limit of the recent SPATS configuration (70 
mPa threshold, > 5 hits per events). The gray band (50-100 mPa 
threshold) around this limit considers uncertainties in absolute noise 
0] and attenuation length Q|. 

192 other experiments. Additionally we show in Fig |5] the 

193 sensitivity for a SPATS-28 detector configuration which 

194 potentially would operate at a lower threshold of 30 mPa 

195 using a new detector readout. 
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